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Abstract 

An analysis of the optical response of a triangular-shaped photonic band-gap prism is presented. 
Numerical simulations have been performed in the framework of multiple-scattering theory, which 
is applied considering spot illumination to avoid diffraction effects. First of all, refractive properties 
in the frequency range below the first TM band-gap are analyzed and compared with the available 
experimental data. It validates the approach employed and supports the predictions obtained in 
the frequency range above the gap. At these high frequencies we found an unusual superprism 
effect characterized by an angle- and frequency-sensitivity of the intensity of outgoing beams. We 
report several representative examples that could be used in device applications. The results are 
interpreted in terms of the corresponding semi-infinite photonic crystal, through the analysis of 
the coupling between external radiation and bulk eigenmodes, using the 2D Layer-Korringa-Kohn- 
Rostoker method. The procedure presented here constitutes a simple but functional alternative to 
the methods used until now with the same purpose. 

PACS numbers: 42.70.Qs, 42.79.Bh, 42.25.Gy, 41.20.Jb 



I. INTRODUCTION 



Since the seminal works of Yablonovitchi and John^, photonic crystals (PCs) have gen- 
erated a great interest, both in fundamental and applied physics. These materials could 
have an important role in the development of the photonic information technology. This 
is especially true in the case of two-dimensional (2D) PCs due to their easy fabrication by 
standard lithographic techniques. On the other hand, recent advances in growing techniques 
are also allowing the fabrication of planar devices based on three-dimensional (3D) systems 3 -. 

The study of optical properties of PC using numerical tool a 4 ' 5 ! 6 ' 7 is an essential procedure 
in order to further the understanding of the physical properties of PCs. The potential 
applications of PCs has motivated the development of efficient theoretical methods that 
allow us to obtain realistic comparisons between theory and experiments. One of the most 
versatile methods is the Multiple Scattering Method (MSM). Its formalism is based on the 
expansion of the 2D (3D) radiation field in cylindrical (spherical) harmonics centered at 
each scattering object, and it has been successfully applied in the analysis of interaction of 
radiation with dielectric and metallic finite 2D and 3D structure o 8 ! 9 ! 10 ! 11 ! 12 ! 13 ! 14 ! 15 ! 16 ' 17 . 

Several authors have analyzed the optical response of structures based on 2D- 
PC o 18 i 19 i 20 i 21 i 22 . Recently, 2D-PCs with non-parallel boundaries have been investigated 
both numerically, by using the finite-difference-time-domain (FDTD) method—, and 
experimentally 2 ^. However, a complete theoretical analysis including the study of the cou- 
pling between incident radiation and bulk eigenmodes is still lacking. 

The goal of this work is to present a comprehensive analysis of the optical properties of 
a photonic bandgap prism by using the Multiple Scattering formalism and other numerical 
techniques. The prism, which has a triangular shape, consists of a hexagonal array of non- 
overlapping dielectric cylinders of infinite height. The analysis is made in the low frequency 
region (below the first TM band-gap), where a comparison with available experimental 
data^ is presented, and also in the frequency domain above the first TM band-gap, where 
anomalous refractive properties are predicted. The results of the MSM simulations are 
interpreted in terms of the corresponding semi-infinite PC. Moreover, in order to get a 
physical insight of the optical response in this structure, we also have made quantitative 
predictions of the coupling strength between the incident light and the bulk eigenmodes 
that can be excited; i.e., the eigenmodes that satisfy the kinetic matching conditions at 
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the air/PC interface. A complete theoretical analysis by means of a combination of plane 
wave expansion (PWE) and 2D-Layer-Korringa-Kohn-Rostoker (2D-LKKR) method, has 
allows us to get accurate predictions which could be confirmed experimentally The good 
agreement found between theoretical results and available experiments supports the different 
approaches employed. 

In our studies we have found an unusual superprism phenomenon, in which the intensity 
of outgoing beams shows a dependence on both the frequency and incident angle of the 
external beam. This differs from the superprism effect reported until now2&, where large 
dependence is found only in outgoing beams directions. Also, it must be emphasized that 
this work provides simple but working guidelines to design functionalities of super-prism 
structures based on a 2D photonic band-gap (PBG) systems. On the other hand, this work 
also shows that the MSM is competitive with the time-demanding FDTD method in order 
to simulate actual 2D structures. 

The paper is organized as follows. In Section II we describe briefly the basis of the MSM 
employed to study the optical response of the prism. The results of numerical simulations 
and the discussion in terms of the semi-infinite model is presented in Section III. Finally, in 
Section IV we summarize the work. 

II. NUMERICAL METHOD 

Here, the MSM is applied to study a finite lattice of dielectric cylinders (of dielectric 
constant e r ) infinitely height embedded in a background of different dielectric constant (e&). 
For a complete description of the numerical algorithm employed see, for example, the paper 
by Ochiai and Sanchez-Dehesa^ and references therein. In brief, the numerical simulations 
are based on the solution of the following self-consistent equatio^ ^ 1 ^ 11 ! 12 ' 1 ^ 14 ' 1 ^ 1 ^ 1 ? : 

C^iatr+E^), (1) 

where the field i/j = E Z (H Z ) for TM (TE) polarization. V? nd in Eq.JD is a column vector 
that contains the components of the multipole expansion of the radiation field induced by 
the cylinder a centered at r a , ip^ represents the external wave, t a is the t-matrix, and G a p 
is the propagator from cylinder /3 to a. 
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Most of the previous simulations using the 2D-MSM employed an incident wave having a 
plane wavefront as external source . However, in the case under study here, a structure with 
non-parallel boundaries, it is necessary to use spot illumination in order to avoid diffrac- 
tion effects. Moreover, spot illumination is essential to make comparisons with available 
experiments^, which use collimated beams in the microwave regime. So far, to the best of 
our knowledge, there have been two different procedures to simulate a collimated beam: (f ) 
the first one defines a finite aperture in front of the sampl e - 10 ' 11 . (2) the second one simulates 
the spot illumination by using a Gaussian bea m 20 i 27 i 28 . We have demonstrated that, for our 
system, a Gaussian beam is more appropriate because it reproduces the measurements fairly 
well. 

The incident Gaussian beam of frequency uj, focusing at r , can be expressed as the 
following very directional solution of the homogeneous Helmholtz equation 29 : 

^ cxt (r) = [ ° dkye- ^^-^, (2) 

J-ko 

where r is the position vector, which in polar coordinates r = (r, 6) . k is the wave vector of 
a plane wave in the beam, k = (k x , k y ), where k x = \Jk,Q — ky, and k^ = uj^ft^jc. Parameter 
D determines the width of the beam. 

Therefore, the £-th component of the external field is given as: 

$g = i% [ \ dO cos fle^W+toW, (3) 

where 



gi (0) = -D jfcg sin 2 fl (4) 
92(0) = k {r a cos {0 + inc - a ) - r cos {0 + inc - O )} - £ (5) 

and (r a , a ) and (ro, Oq) are the polar coordinates of the a-th cylinder and the focusing point, 
respectively. 

We are interested in the far field response of the structure. It can be easily found that in 
the far field regime^: 

^(r) -> ( a ° {9) etk0r/V¥ if - f " °- f (6) 

a o (0 — n) e^ lk ° r / y/r otherwise, 
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where a (0) = y/2nk e~ 1 ^ e~ Dk ° sm2 9 cos9. 

The magnitude of interest in this problem is the intensity at infinity (J s ). It is defined 
by the Poynting vector flux per unit angle calculated at far field: 

H W») + /.W if -!<«<! 

l/.(«)| 2 - M»)| 2 otherwise, 
where f s {0) is the scattering amplitude, which is obtained from Eq.Q as a linear combination 
afVffl (Ref. 15). 



III. RESULTS AND DISCUSSION 



The optical response to TM-polarized light of a 2D photonic band-gap prism based on 
a hexagonal array of 1540 alumina rods (e r =8.9) in air (e?, =1.0) is studied. In order to 
compare with the experiments performed by Lin et alM the radius of cylinders is set to 
r =0.188 a, where a is the lattice parameter. The external shape of the prism is triangular: 
the three boundaries of the prism are composed of 55 cylinders each one, and they are 
parallel to the TK direction. A scheme of the structure and the definitions of incident angle 
{Oinc) and deviation angle (S) are shown in Fig. 1. The system of reference used along the 
calculations is also plotted. 

The photonic band structure for TM-polarized modes of the corresponding infinite lattice, 
which has been calculated by means of a plane wave expansion 31 , is plotted in Fig. 2, where 
the frequencies are given in reduced units. A factor of 2nc/a (c is the light velocity in air) 
must be used to obtain the absolute values. So, in the rest of the paper the frequency is 
given in reduced units. Results in Fig. 2 indicate that a complete band-gap exists between 
u L =0.325 and u v =0.497. 

First of all, frequencies inside the band-gap were examined in order to check whether the 
bulk properties of the hexagonal lattice were already achieved. Let us consider a frequency 
slightly greater than ul, for instance, 0.33. The electric pattern computed for the case of 
normal incidence (0j nc =0°), D =20 a 2 , and ro=0 is visualized in Fig. 3. It can be seen that 
no light propagates inside the prism and, therefore, this frequency is inside the photonic 
band-gap, which is in agreement with the photonic band structure shown in Fig. 2. 

In what follows we report the optical response in two different regimes, below and above 
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the first photonic TM band-gap. We obtained converged results by considering I = 1 as 
the maximum angular momentum for frequencies below the band-gap, and i = 2 for the 
frequencies considered above the band-gap. 

A. Behavior at frequencies below the first TM band-gap 

In the very low frequency region the equifrequencial curves are circles and, therefore, it 
is possible to define an effective refraction index independent of incident angle^. 

It can be stated that, by using the MSM and for frequencies below the first TM band-gap, 
it is possible to accurately find the deviation angle of the beam transmitted across the prism 
only by assuming certain ro j^O and D values large enough^. 

In addition, the asymmetry observed in the transmitted peaks can be explained by the 
Fresnel factor because of the fact that the incident beam is represented by superposition of 
plane waves with different incident angles and different weights [Eq.(2)]. 

A procedure to get the effective refractive index (n e ff) from the MSM is to fit the numer- 
ical results (9 inc , S) to the expression obtained applying Snell's law in the case of non-parallel 
boundaries structure with a homogeneous n e //: 

5 = Oinc+ sin" 1 ((sina)(n 2 eff - sin 2 9 inc )* - sin 9 inc cos a) - a, (8) 

where a is the prism angle, which is 60° for the triangular prism under analysis. Figure 
4 shows the dependence of 5 on 9 inc for several angles at u;=0.26 (black dots). The solid 
line corresponds to a fit of MSM data to Eq.(JHJ). From this fit n e // =1.54 is obtained. 
Experimentally, Lin et a/.— reported n e ff «1.57 (dashed line in Fig. 4). Therefore, our 
simulations based on the MSM confirms that at this low frequency (u=0.26) the refractive 
index does not depends on 9i nc and, moreover, its value fairly agrees with experimental data 
(about 2% of relative error). 

The previous considerations allow us to analyze theoretically the optical response for 
all frequencies below the first band gap. Unfortunately, we found that the contribution of 
the transmitted beam for frequencies near the band edge are not clearly observed in the 
calculations (its transmission is almost zero). In other words, the scattered intensity is 
mainly determined by the shape of the prism. The result is summarized in Fig. 5(a), where 
n e ff obtained from the MSM simulations (black dots) are compared with experimental data 
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(white circles) and with a simplified model (solid line) that considers light propagation across 
single air/PC interfaces (see Appendix A). The last approach predicts a sharp decreasing 
of n e ff near the band edge, which, unfortunately, is confirmed neither by the MSM nor by 
the experimental set up due to their intrinsic limitations. Figure 5(b) shows the behavior of 
n e ff as a function of 6 inc . It can be observed that, as it is expected, n e ff does not depend on 
9 inc at low frequencies, but near the band edge, the critical frequency at which n e // starts 
to decrease depends on 9 inc . This phenomenon is a result of the shape of the equifrequencial 
curves near the band edges. This prediction should stimulate more accurate experimental 
works looking for its demonstration. 

Finally, let us remark that small discrepancies observed between experimental and numer- 
ical results can be mainly attributed to the fact that the actual beam used experimentally 
is not well enough simulated by the 2D numerical beam introduced in the simulations by 
the MSM. In other words, the spot size and collimation of a 2D Gaussian beam cannot 
be controlled to reproduce accurately the properties of the actual beam. This is the main 
drawback of the simulations based on the MSM, which produces fictitious diffraction effects. 
Fortunately, for higher frequencies these effects disappear, and it is possible to characterize 
interesting phenomena as it is reported in what follows. 

B. Behavior above the first TM band-gap 

In order to study the optical response in frequencies above the first complete TM band- 
gap, we have started by analyzing the case of infinite structures and have chosen the ranges 
of frequency and incident angle at which anomalous phenomena are expected when light 
crosses an air/PC interface. Equifrequencial curves between 0.509 and 0.627 are plotted in 
Fig. 6(a). Notice that near the curves associated to frequency 0.537 [thick lines in Fig. 
6(a)] the equifrequencial curves consist of two separated regions: a hexagonal shaped region 
(centered at T point) and a triangular shaped region (centered at K point). Therefore, it is 
expected that such anisotropic equifrequencial regions lead to anomalous response of a finite 
structure based on the same lattice periodicity. 

The procedure described in Appendix A allows us to compute the propagating angle 
(8 pro ) inside the PC as a function of the incident angle for a beam crossing the air/PC 
interface at u> =0.537. The result of this calculation, which is depicted in Fig. 6(b), shows 
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that above a critical angle (9 C ) the incident beam splits into two. They correspond to 
the curves of frequency 0.537 which have group velocity with positive y— component. The 
great dependence of 9 pro on 9 inc of these solutions from 20° to 25° degrees and from 25° to 
40° is associated to the rounded corners of the hexagonal- and triangular-shaped regions, 
respectively. 

A comprehensive analysis of the scattering properties in the frequency range between 0.5 
and 0.6 and for incident angles between 20° and 40° has been performed by using the MSM. 
This analysis has allowed us to select the most interesting phenomena. 

Figure 7(a) shows the intensity map as a function of S and 9i nc computed with cj=0.537. 
The parameters of the Gaussian beam employed in both simulations^ are D =30a 2 and 
y =10a . As can be seen from this figure, we do not obtain a great dependence of 5 as a 
function of 6i nc as could be expected in the usual definition of superprism effect. Instead, 
we obtain an angular sensitivity of the transmitted beam intensities. In order to prove 
clearly this anomalous effect, Fig. 7(b) shows the comparison between intensity spectrum 
corresponding to the incident angle with largest transmission (6^=22.8°), and the spectra 
corresponding to incident angles that differ in 3.2°. As it is shown in this figure, a small 
change in 9 inc implies a dramatic fall in the transmittance. 

In addition, it is interesting to know the dependence of I s on uj and 5 for #j nc =22.8°. 
The result of this calculation is shown in Fig. 8(a). This second transmitted beam only 
appears in a very narrow range of frequencies; i.e., from 0.533 up to 0.582 approximately. In 
order to show quantitatively this effect, Figure 8(b) shows a sharp change in the intensity 
features of the spectra for very small variations (about 0.7 %) of frequency. Again we obtain 
a novel refractive effect characterized by the large frequency-sensitivity in the intensity of 
the outgoing beams. 

The phenomena mentioned previously suggest the above configuration as a sensitive beam 
splitter, both in angle- and frequency regime. 

A similar analysis to the one described previously was carried out for another frequency 
of interest, c<j=0.539. The corresponding intensity map, with 9i nc and 5 as variables, is 
plotted in Fig. 9(a). From this figure we have selected the angle #j nc =40° as the most 
interesting case. Figure 9(b) shows its corresponding intensity spectrum compared with the 
spectra associated to two angles differing in 5.2°. As can be seen from this figure, the main 
transmitted beam found at #j r j C =40° disappears completely for #j nc =34.8°, while it is almost 
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the same for #j nc =45.2°. 

Also we have made a frequency analysis by fixing 0j nc at 40°. The results are summarized 
in Fig. 10(a), where it is seen that the main transmitted beam only exists in two separated 
frequency regions. Particularly, Fig. 10(b) shows the behavior in the region centered at 
u;=0.539. Notice that a change in frequency of about 4% implies a sharp decrease in the 
transmission of the incident radiation through the prism. Therefore, this configuration could 
be suitable for the design of a PC filter for a short range of frequencies. Similar conclusion 
could be obtained for frequencies above 0.58, where the prism is almost transparent and the 
transmitted beam undergoes a small deviation angle (5 ~ —20°). 

Physical insight of the mechanism originating the phenomena observed in Figs. 7, 8, 9, 
and 10, can be obtained from the 2D-LKKR metho d 35 ^ 36 . In fact, we have implemented a 
2D-version of the procedure proposed by two of us in a previous work2i (see Appendix B). 
In this way, it is possible to obtain the branching ratio (BR) of the different bulk eigenstates 
of the PC that can couple to the incident light. 

In the following, we will refer to the configurations defined by (u, 6 , inc ) = (0.537, 22.8°) 
and (u, inc ) = (O.539, 40.0°) 1 and case 2, respectively. 

In case 1 (see Fig. 11), simulations using the MSM predict that in addition to the 
reflected beam (R) at the front air/PC interface (133.27°), two transmitted beams appear: 
one of them at -105.82° (Tj) and another at 14.42°(T 2 ). In this case, if we apply the kinetic 
matching conditions to the frontal air/PC interface, we obtain that only one eigenstate with 
group velocity inward the PC couples to the incident light. The perpendicular component, 
k±, of this eigenstate wave vector (k pro ) to the boundary is -0.3592 (in units of 2n/a) [see 
inset in Fig. 11(b)]. A similar analysis can be performed at boundary 2 in Fig. 11(a). Since 
boundaries 1 and 2 are not parallel, it is necessary to make a 60° rotation of k pro before 
applying the matching conditions. As a result of this analysis, we conclude that an output 
beam should appear at 5 =-105.48°, in good agreement with the MSM-based simulations. 
Also we have to consider the reflection in boundary 2 and the following refraction in boundary 
3. For this path the predicted deviation angle is 14.4°, and we expect small intensity for this 
peak since it comes from considering two refractions and one reflection. On the other hand, 
from the 2D-LKKR it is possible to obtain the reflectance of the corresponding semi-infinite 
PC at this incident angle. The reflectance is about 0.544. This result is supported by the 
MSM, where the intensity of the peak associated with the reflected beam (R in Fig. 11) is 
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related to the intensity of the incident beam by that factor. 

The field pattern and transmission spectrum obtained from the MSM-based simulations 
for case 2 (uo =0.539 and 9 inc =40°) are plotted in Figs. 12(a) and 12(b), respectively. A very 
small reflected peak (R) appears at 101.16°. The main transmitted beam (Ti) with negative 
refraction index appears at -146.16°, while two secondary transmitted beams are seen at - 
93.24 (T 2 ) and -19.98° (T 3 ). This result can be understood in similar terms to the preceding 
case. In other words, by analyzing the light transmission through single interfaces. Now 
we have two eigenstates satisfying the kinetic matching conditions at the frontal interface; 
their k± values inside the PC, in units of 2%/a, are 0.5095 and -0.1237 [see inset in Fig. 
12(b)], and the corresponding BRs are 0.923 and 0.077, respectively. The reflectance of the 
corresponding semi-infinite PC is about 0.085. This value of the reflectance is in accordance 
with the small intensity peak of the reflected beam in comparison with the main transmitted 
beam. The matching conditions applied to boundary 2 in Fig. 12(b) concludes that the first 
outgoing beam should appear at 5 =-145.75°, in good agreement with the MSM results. If 
we also take into account the beam with negative k±, we obtain that a second outgoing beam 
appears at 5 =-93.0° (also in agreement with Fig. 12). The great difference between the two 
intensities of these two refracted beams could be estimated using the BR values. Moreover, 
for the path corresponding to the light reflected at the boundary 2 and refracted at boundary 
3, the predicted deviation angle is -20.0° and also we expect a secondary refracted beam, 
since the dominant eigenstate in boundary 2 is the refraction. 

Table I summarizes the results regarding the directionality of the scattered light. It 
demonstrates that the deviation angle of the transmitted light through non-parallel bound- 
aries can be predicted by a simplified model based on the light transmission across semi- 
infinite air/PC interfaces. Its results agree fairly well with the complete simulation of the 
finite system performed by MSM-based simulations. In other words, the combination of 
the Hellmann-Feynmann theorem in the plane wave expansion together with the 2D-LKKR 
method is enough to analyze the optical response of this kind of structures in the high 
frequency region. 
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TABLE I: Comparison of the deviation angles (in degrees) of beams transmitted across the prism. 
The results under the columns 'Model' are obtained by the procedure described in Appendix A. 



lj=0.537, e inc =22.8° ^=0.539, 9 inc =40° 

MSM Model MSM Model 

Ti beam -105.82 -105.48 -146.16 -145.75 

T 2 beam 14.42 14.40 -93.28 -93.00 

T 3 beam -19.98 -20.00 



IV. SUMMARY 



In this work we have analyzed the optical response of a triangular prism made of a 
hexagonal array of dielectric cylinders in air. The approach employed, based on the Multiple 
Scattering formalism, has been supported by the available experimental data existing at 
frequencies below the first TM band-gap. Also, we have predicted anomalous superprism 
phenomena above the first TM band-gap. A comprehensive analysis in this high frequency 
region has been performed by combining the PWE method and the 2D-LKKR method. We 
have demonstrated that, if the prism is large enough (diffraction effects are negligible), a 
procedure based in the analysis of the corresponding semi-infinite system agrees very well 
with the simulation of the complete finite system by the MSM. We have described two cases 
with possible device application. The results presented in this work support the applied 
procedure as an efficient numerical tool in the understanding of the optical response in 
2D-PC structures with non-parallel boundaries. We hope that this work motivates further 
experiments that confirm its predictions. Finally, the procedure here employed has been also 
extended to study analogously 3D prisms composed of PCs having a strong anisotropy in its 
photonic dispersion relation, and the results will be presented in a forthcoming publication. 
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APPENDIX A: LIGHT TRANSMISSION ACROSS A SINGLE AIR/PC INTER- 
FACE 

Let us assume a semi-infinite 2D PC defined in the region y >0, and air in the another 
semi-plane. In order to know which eigenstates of the PC can be excited by an external 
light, it is important to project the bulk photonic band structure on the surface Brillouin 
zone (SBZ). A typical projected band diagram is shown in Fig. 13. This figure represents 
the band structure projected on the YK boundary surface, which is the common boundary 
surface of the triangular prism under study. In Fig. 13 the shaded regions define the 
allowed bulk eigenstates and the thick black line defines the light line. This line defines 
the critical frequency above which external plane waves are allowed to exist. The white 
regions represent pseudogaps in the dispersion relation. Therefore, external beams having 
frequencies in a pseudogap cannot enter into the PC, as it is shown in Fig. 3. 

Another criterion of coupling arises from the group velocity of the bulk eigenmodes. Only 
the bulk eigenstates with positive v y can couple with the incident light. Those eigenstates 
with v y < should be physically neglected. 

In summary, the kinetic matching condition between the bulk eigenmodes and the incident 
light are the following: (1) The frequency and the reduced wave vector in the SBZ are 
conserved, (2) The bulk eigenmodes are above the light line, and (3) v y is positive. Once 
these conditions are satisfied, the incident angles 6 inc and the propagating angles 9 pro , which 
are given by 

k = — (sin 9 inc , cos6 inc ), (Al) 
c 

v = \v |(sin# pro , cos9 pro ), (A2) 
can be related through an effective refraction index (n e ff) defined as: 

sin(6 inc ) = n e //sin ^arctan (^) 

where (v x ,v y ) are the components of the group velocity, which are computed directly from 
the dispersion relations by using the Hellmann-Feynman theorem. 
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APPENDIX B: BRANCHING RATIO 



Figure 2 shows that at high frequencies several eigenmodes can be kinetically matched to 
the external light. Hence, it is necessary to know quantitatively what the coupling strength 
is between the incident light and the bulk eigenstates of the 2D PC. Here, the 2D layer-KKR 
method is used for that purpose in a similar manner as was employed its 3D counterpart in 
the analysis of the superprism effect in opals^i. 

We have considered a slab large enough for the deepest region in the slab to be regarded 
as a superposition of the actual bulk eigenstates a; i.e., the field at this center region is 
assumed not be affected by the finite size of the slab. In the 2D-LKKR method three 
different regions can be clearly separated: (1) Left region, defined by the scattering matrix 
Ql- (2) Center or deepest region. (3) And right region, whose scattering matrix is Qr. If 
we consider that a plane wave impinges the slab from the left, the following relationships 
are satisfied 34 ' 35 : 



E 



\ E re/ 



E 

e: 



tr 



)- 


(q l ++ 




( 




{Q L - 


Q'+j 
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( Q R ++ 


Ql-\ 


( 




K Q R - 








o 



(Bl) 



(B2) 



where Ej nc , E re /, E tr are column vectors which contain the Fourier components of the 
incident, reflected, transmitted fields, respectively. E^ corresponds to the field in the inner 
void region of the slab, '+' denotes propagation from left to right and '— ' from right to left. 
From Eq. (|Bl|) and (|B2|) it is easy to obtain: 



E+ = (1 - Q L + _Q R + Y l Q L ++ E inc (B3) 
E^ = Q R + (1 - Q^Q^+r 1 Q L ++ E inc (B4) 

Let us consider the a-th left eigenstate (V ai ) of the transfer matrix which relates the 
Fourier components of the electric between two consecutive voids in a infinite PC. It is 
possible to define the projection of the eigenstate in the deepest void of the slab in the true 
a-th eigenstate through the coefficients c a defined by: 
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c a = (V+ L )' E+ + (V- L )* EJ (B5) 
These coefficients allows to define the BR corresponding to the bulk eigenstates: 







2 




c a >\ 


2 



where the summation is over all the bulk states satisfying the kinetic matching condition. So, 
BR Q gives an estimation of the coupling of the external light with the cc-th bulk eigenstates 
in the 2D PC, and it is normalized to the total coupling of eigenstates available for light 
transmission. 
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FIG. 1: Schematic view of the triangular prism under study and parameter definitions. The 
reference system used in the calculation is also plotted. 

FIG. 2: Photonic band structure for TM-polarized modes corresponding to a hexagonal lattice of 
dielectric cylinders (e r =8.9) in air (ef,=1.0). The cylinder radius is r=0.188a, where a is the lattice 
constant. The shaded region defines the complete band gap. 

FIG. 3: Electric field pattern, E z (x,y) generated by a Gaussian beam, with the parameters 
D=20a 2 , ro=0 and w=0.33, impinging on the prism. 

FIG. 4: Dependence of deviation angle (<5) on incident angle (0j nc ) for uj =0.26 (solid points). Solid 
line is obtained by fitting Eq.(jHJ. n e //=1.54 is obtained from this fit. Dashed line corresponds to 
experimental results. 

FIG. 5: (a) Behavior of the effective refractive index (n e ff) as a function of frequency. Hollow 
points corresponds to experimental data, solid points to MSM simulations and solid line to the 
simplified model (see text), (b) Dependence of n e ff as a function of the frequency for different 
incident angles (0i nc ) computed by the simplified model. 

FIG. 6: (a) Equifrequencial curves for the photonic TM-bands in Fig. 2. Curves for frequencies 
between 0.509 and 0.627 are represented. Solid line corresponds to w=0.537. (b) Dependence of 
the propagation angle on the incident angle considering u =0.537. Dashed line is the result of 
Snell's law calculated by considering an uniform medium with a spatially averaged refractive index 

Wave — 1.419. 

FIG. 7: (a) Intensity map as a function of the deviation angle ((5) and the incident angle (^mc) 
computed with a;=0.537. White corresponds to zero values and black to the maximum intensity. 
Dashed line defines the incident angle in which a maximum of the transmittance is found, (b) 
Comparison between the intensity spectra corresponding to the same frequency (a;=0.537) and 
three different angles. 'R' labels the peak corresponding to the reflected intensity. Notice the 
angular-sensitivity of intensity features. 
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FIG. 8: (a) Intensity map as a function of the deviation angle (5) and the frequency (a;) computed 
with 9i nc =22.8° . White corresponds to zero values and black to the maximum intensity. Dashed 
line shows the frequency in which a maximum of the transmittance is found, (b) Comparison 
between the intensity spectra corresponding to the same incident angle (6*j nc =22.8 ) and three 
different frequencies. 'R' labels the peak corresponding to the reflected intensity. Notice the 
frequency-sensitivity of intensity features. 

FIG. 9: (a) Intensity map as a function of the deviation angle (<5) and the incident angle (^nc) 
computed with a;=0.539. White corresponds to zero values and black to the maximum intensity, 
(b) Comparison between the intensity spectra corresponding to the same frequency (u;=0.539) and 
three different angles. 'R' labels the peak corresponding to the reflected intensity. Notice the 
angular-sensitivity of intensity features between the top and the center spectra. 

FIG. 10: (a) Intensity map as a function of the deviation angle (S) and the frequency (oo) computed 
with #j nc =40.0°. White corresponds to zero values and black to the maximum intensity. Dashed 
line defines the frequency in which a maximum of the transmittance is found, (b) Comparison 
between the intensity spectra corresponding to the same incident angle (#j nc =40.0°) and three 
different frequencies. Notice the great frequency-sensitivity of intensity features. 

FIG. 11: (a) Electric field pattern and (b) transmission spectra generated from a Gaussian beam 
with lo =0.537 and #j nc =22.8° . White arrows indicate the direction of the group velocity computed 
by the semi- infinite model (see Appendix A). The solutions to k± and the BR value are shown in 
the inset. 

FIG. 12: (a) Electric field pattern and (b) transmission spectra generated from a Gaussian beam 
with oj =0.539 and 9i nc =40°. White arrows indicate the direction of the group velocity computed 
by the semi- infinite model (see Appendix A). The solutions to k± and the corresponding BRs values 
are shown in the inset. 

FIG. 13: Projection of the photonic band structure shown in Fig. 2 on the surface Brillouin zone of 
the TK surface. The shaded regions defines the bulk eigenmodes. The frequencies of the external 
radiation field are defined above the light line (black line) , u > ck\ | . 



17 




Fig. 2. J. Bravo-Abad et al. 




Fig. 3. J. Bravo-Abad et al. 




Fig. 4. J. Bravo-Abad et al 



1.8-r 

1.7- 
1.6- 
1.5- 
1.4- 
1.3- 
1.2- 



(a) 



PB( g 



o Lin et al. 
• MSM 

Semi-infinite 



o cm 



0.18 0.20 0.22 0.24 0.26 0.28 0.30 0.32 0.34 

Frequency (red. units) 




0.18 0.20 0.22 0.24 0.26 0.28 0.30 0.32 0.34 

Frequency (red. units) 



Fig. 5. J. Bravo-Abad et al. 



(a) 





0.6- 




0.4- 




0.2- 


■ a 




C\J 


0.0- 












-0.2- 




-0.4- 




-0.6- 




0.6400 
0.6200 
0.6000 
0.5800 
0.5600 
0.5400 
0.5370 
0.5200 
0.5000 
0.4800 



-0.6 -0.4 -0.2 0.0 0.2 0.4 0.6 

k a/2% 

X 



CD 
CD 



I "20 
-40 
-60 
-80 



(b) 



80 
60 
40 
20 



CD o -[•^••••••••• ## 



T ' 1 ' 1 ' 1 ' 1 ' 

• • Semi-infinite model 
• Snell's law 



9 



10 20 30 40 

e inc( de 9-) 



50 



60 



Fig. 6. J. Bravo-Abad et al. 



(a) 



150! 
100- 
50- 



C Oh 

to 

-50 k 



-100- 
-150- 



20 24 28 32 36 40 

0. (°) 

inc x ' 

(b) 



15 
10 
5 

£ 
c 15 

I 10 

_» 
15 

10 

5 





j i_ 



0. =19.6° 

inc 



9. =22.8° 



il i L 



1 



' I V i L 



9. =26.0° 



I i L 



R 



J i L 



R 



R i 



' A- . | , L 



150 -100 -50 50 100 150 

8 



Fig. 7. J. Bravo- Abad et al. 



(a) 



o 

CO 



150 : 

100- 
50- 
0- 
-50- 
100- 
150 - 

0.50 0.52 0.54 0.56 0.58 0.60 
Frequency (red. units) 



B 



3 

co 



15 







(b) 



II I 1 I 1 I 1 I 1 I 1 I 1 

' co=0.533 R 




1 


| i | i | i | i | i | , 

co=0.537 p 

1 i A 


■ i 1 i 1 i 1 i 1 i 1 i 1 
: co=0.541 pj 

- \ 1 1 ^ 1 1 1 1 1 1 1 


1 



150 -100 -50 50 100 150 
5 (deg.) 



Fig. 8. J. Bravo-Abad et al. 



(a) 






15- 




10- 




5- 




0- 




15- 


"E 


10- 


=3 


5- 


_Q 


0- 


15- 




10- 
5- 
0^ 







6. =34.8° 

inc 


R 




1 i A A 1 


, 1„ 






e. =40.0° 

inc 










R 






i A i *. i 








6 =45.2° 

inc 




u 






R 










-150 -100 


-50 50 


100 150 






5(°) 





Fig. 9. J. Bravo-Abad et al. 



o 

CO 



-100 - 



15- 
10- 
5- 



oT 
■E 15h 
3 10 



5h 


15- 

10- 

5- 





(a) 




0.50 0.52 0.54 0.56 0.58 0.60 
Frequency (red. units) 



(b) 



co=0.519 



co=0.539 



, i A. 



J C^-, L 



co=0.559 



J L-a L 



150 -100 -50 50 100 150 

6(°) 



Fig. 10. J. Bravo-Abad et al. 



(a) 




-150 -100 -50 50 100 150 



8(°) 

Fig. 1 1 . J. Bravo- Abad et al. 



40- 



20- 



0- 



H — > 

-Q 



-20- 



-40- 




' O O c> o» «o O-c'-O 
© o - c - <■ © * 0» C> • o o 

=- *<> ^<S>/^ C^- O o o • C3 
> O • o « •>-> •o««c?«o © • 

O o © O o COO 

■=>•<> i<>>pjcy^^ ° • • 

.ot<>io <>. ^7 : * 



T 



25 
20 
15 
10 
5 




-40 -20 

(b) 



20 40 



T 



0.550 



€ 0.545 



1 0.540 



0.535 



0.530 





7.7% 




92.3% 










1 * 
















ri 


J \ \ 




-0.6 


-0.4 -0.2 


0.0 


0.2 0.4 


0.6 



T 



R 



150 -100 -50 50 100 150 

5 (deg.) 



Fig. 12. J. Bravo-Abad et al. 




Wa/2% 



Fig. 13. J. Bravo-Abad et al. 



